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ABSTRACT 

In paper I, we showed that time-dependent general relativistic magnetohydrody- 
namic (GRMHD) numerical models of accretion disks, although being highly turbu- 
lent, have surprisingly simple electromagnetic properties. In particular, the toroidal 
current density in the disk takes the form dl^/dr oc r -5 / 4 . Guided by this simplicity, 
we use a time-dependent general relativistic force- free electrodynamics (GRFFE) code 
to study an idealized problem in which the accretion disk is replaced by an infinitely 
thin rotating equatorial current sheet. We consider both an r -5 / 4 current profile and 
an r _1 profile, the latter corresponding to the paraboloidal model of Blandford & 
Znajek (1977). The force-free magnetosphere we obtain with the r -5 / 4 current sheet 
matches remarkably well to the Poynting-dominated jet seen in GRMHD numerical 
models. By comparing to the non-rotating force-free model studied in paper I, rotation 
is seen to lead to mild decollimation of the jet suggesting that hoop-stress forces nearly 
cancel centrifugal forces. In order to study the process that generates the corona and 
disk wind and destroys the ordered field in the corona in GRMHD numerical models, 
the force-free field with the r -5 / 4 current distribution is embedded in an accretion 
disk and followed in a GRMHD simulation. Field at high latitudes is continuously 
transported to larger radii leaving a corona with only disordered field, while in the 
equator the turbulent field is accreted. Reconnection and magnetic stresses contribute 
to a magnetized, thermal wind without the aid of an ordered field threading the disk. 

Key words: accretion disks, black hole physics, galaxies: jets, gamma rays: bursts, 
X-rays : bursts 



1 INTRODUCTION 

Accretion disks around black holes have been re- 
cently studied using general relativistic magnetohydro- 
dynamic (GRMHD) simulations that have been used 
to determine the structure of the disk, corona, wind, 
and Poynting-dominated jet Jibe Villiers. Hawlev. fc Krolikl 
120031 ; IMcKinnev fc Gammiel [20041 ') . This paper continues a 
study of the coupling between these components of the ac- 
cretion flow. 

The simulated accretion flows develop vigorous non- 
linear tu rbulence driven by the magn eto-rotational in- 
stability (jBalbus fc Hawlevl Tl99ll . Il998l h One might sus- 
pect that all fluid and electromagnetic quantities in the 
dis k would therefore be chaoti c. However, as we showed 
in IMcKinnev fc Naravanl l|2006j ). the vertically integrated 
toroidal current exhibits a smooth and simple behavior 
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ijMcKinnev fc N aravan 2006). Specifically, the toroidal cur- 
rent 1$ enclosed inside radius r closely follows a power-law 
profile, 



dr r- 



(1) 



with v = 3/4. This scaling is found to hold in a time- 
averaged sense but also at each instant of time. It is also 
independent of the black hole spin or the initial cond itions 
used for the simulation. IMcKinnev fc Naravanl (|2006l ) ar- 
gued that the presence of the simple current distribution is 
the reason why GRMHD models have simple, nearly steady, 
collimated, magnetically-dominated polar regions which are 
qualitatively similar in structure to t he force-free jet in the 
Blandford-Znajek (1977, BZ) model iMcKinnev fc Gammiel 
l2004l : lHawlev fc Krolikl 120061 ; lMcKinnevll2006d\ ~ 

The fact that the current distribution in GRMHD 
models is simple suggests that the accretion disk might, 
at some level of approximation, be idealised as a simple 
boundary condition in the equatorial plane. Motivated by 
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this possibility, we consider general relativistic force-free 
electrodynamical (GRFFE) models with equatorial bound- 
ary conditions derived from GRMHD numerical models 
of turbulent accretion disks. Treating the accretion disk 
(or stellar surface) as a boundary condition is a com- 
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magnetosphere of the accretion disk leads to simple 
tractable equations that can be solved analytically (or quasi- 
analytically) in certain regimes, such a s for slowly rotating 
black holes (IBlandford fe Znaiek)ll977l ) or for special field 
geometries l|Uzdensky||2004 120051 ) . 

The most famous GRFFE model that treats the disk as 
a boundary condition is th e Blandford-Znajek (BZ) model 
l|Blandford fe Znaiekll 19771 ) , which involves a stationary disk 
with a fixed toroidal current density. BZ presented two solu- 
tions for slowly spinning black holes with an equatorial cur- 
rent sheet. The first solution corresponds to a split-monopole 
field which is centered on the black hole; this solution is a 
general relativistic ext ension o f the force-free pulsar mag- 
netosphere solution of iMichell l| 19731 ) and has a power-law 
disk current of the form given in equation fl} with v — 0. 
T he second solution corresponds to a paraboloidal field, as 
in lBlandfordl (|l976l ). and has v = 1. 

While the monopole solution is not applicable to colli- 
mating jets, BZ's paraboloidal solution is a potentially inter- 
esting model of astrophysical jets. However, the paraboloidal 
solution is incomplete because BZ do not ensure force bal- 
ance across the boundary between field lines threading the 
black hole and those threading the equatorial accretion disk. 
BZ identified this interface between the black hole and ac- 
cretion disk as the transition region. Their estimate of the 
energy output of the paraboloidal solution may be quali- 
tatively incorrect since a significant amount of the energy 
comes from the transition regio n where the disk cont a cts th e 
ergosphere. Since t he work of IBlandford fe Znaiekl l| 19771 ). 
iMacDonaldl l| 19841 ) ; iGhosh fe Abramowiczl (|l997l ) considered 
better-motivated force-free solutions in order to estimate 
the flux ratio of the hole and disk but did not account 
for the transi t ion re gion between the black hole and disk. 
iMeliani et alj (|2006l ) recently developed some small-angle 
approxima tions for the pol ar jet for a Schwarzschild black 
hole, while lUzdenskvl l|2005l ) studied closed field (no jet) ge- 
ometries for a Kerr black hole. 

The difficulty of finding solutions to the GRFFE and 
GRMHD equations of motion has spawned a significant in- 
terest in developing time-dependent codes to solve these 
equations of motion. Of relevance to this paper is that 
time-dependent force-free codes f o und th a t the monopo- 
lar BZ solution i|Komissarovl|200ll, l2002allbl. |2004 and the 
split-monopole BZ solution (|McKinnevl 12006a ) are stable. 
More recently force-free and MHD (in the highly mag- 
netized limit) codes have been used to successfully study 
the pulsar magnetosph ere in the regime that all field lines 
cross the light cylinder (IKomissarov 2006; M cKinnevll2006bl ; 
iBucciantini et al.ll2006l ; IspitkovskvfeoOrj) . 

Our first objective in this paper is to use a time- 



dependent GRFFE code (|McKinnevl l2006al ) to evolve the 
black hole + equatorial current sheet system with disk ro- 
tation and an arbitrarily rapidly rotating black hole to self- 
consistently obtain the complete paraboloidal-type field so- 
lution. 

iMcKinnev fe Naravanl ([2006) found that the current 
distribution and poloidal funnel field geometry in GRMHD 
simulations of accretion disks are not consistent with BZ's 
paraboloidal model with v — 1 or BZ's monopolar model 
with v = 0. Instead, the best fit to the current distribu- 
tion and the poloidal funnel field geometry is obtained with 
v = 3/4. Despite the apparent similarity between the v = 1 
and v = 3/4 solutions, they are qualitatively different. For 
example, the v = 1 solution has a logarithmic divergence 
as a function of radius for the enclosed toroidal current. On 
the other hand, the = 3/4 solution has a finite enclosed 
toroidal current. Given the similarity between the GRMHD 
simulations and the v = 3/4 current distribution and result- 
ing funnel field geometry, our second objective is to solve 
the time-dependent GRFFE equations of motion for this 
GRMHD-motivated model with v = 3/4. 

In order to construct a force-free model that closely 
matches the GRMHD simulations, we require not only the 
current density in the equatorial sheet but also the angular 
frequency of the field lines at the equator. The simpl e force - 
free models we described in Mc Kinney fe Naravanl (|2006T ) 
had no rotation, but here we are interested in modeling the 
effects of rotation. For this purpose we use the results of 
GRMHD simulations to construct a reasonable model of the 
field rotation frequency for th e entire radial range from the 
horizon out well into the disk. I McKinney fe Naravanl (|2006T ) 
found that the disk field angular frequency is Keplerian at 
large radii but makes a smoothly transition near the black 
hole to VL f /VLh ~ 0.5 for a/M > 0.4 and Q F /Q H ~ 1 for 
a/M < 0.4, where Qh is the angular frequency of the hole. A 
GRFFE solution with this dependence for £If and a current 
distribution of v — 3/4 would be the closest a GRFFE model 
could come to the GRMHD simulations without including 
the vertical structure of the disk or the matter itself. Also, 
with £If in the transition region modelled accurately, the 
energy output of the disk and the black hole can be measured 
and compared. We study and compare the jet structure in 
the idealized GRFFE model w ith that found in GRMHD 
simulations of accretion disks (|McKinnevl l2006cl ) and find 
remarkable agreement between the two. 

Our final objective is to demonstrate how the v = 3/4 
force- free solution changes as a result of the presence of mat- 
ter and the magneto-rotational instability. We choose the 
v = 3/4 type stationary field geometry as the initial condi- 
tions and carry out a GRMHD simulation with an accretion 
disk superimposed onto the field. We find that the model de- 
velops a turbulent disordered disk, corona and wind, but the 
Poynting-dominated jet power output is unchanged relative 
to the GRFFE model. 

Paper Outline 

In section [2] we outline the numerical setup and the method 
used to obtain force-free solutions. In section[3] we study the 
v = 1 paraboloidal force- free solution. In section [3] we study 
the v = 3/4 force- free solution using GRMHD-motivated 
boundary conditions for the disk. In section [5] we discuss 



© 2006 RAS, MNRAS 000,1111201 



Disk-Jet Coupling in Black Hole Accretion Systems II: FFE Models 3 



a GRMHD numerical model that is initialized with an ac- 
cretion disk embedded with the v — 3/4 force- free field. We 
follow the destruction of the force-free field associated with 
the corona in GRMHD models. In section [6j we discuss the 
limitations of our calculations. Finally, in section [7] we dis- 
cuss the results and conclude. 



Units and Notation 

The units in this paper have GM = c = 1, which sets the 
scale of length (r g = GM/c 2 ) and time (t g = GM/c 3 ). The 
density scale is arbitrary, so values of the density are nor- 
malized by some fiducial field strength. In order to convert 
to a physical value of the density for a given mass accretion 
rate, one must use the field strength as a function of black 
hole spin and mass accre tion rate as given b y G RMHD mod- 
els such as described in IMcKinnevI l|2005alBB, Ho06c ) . The 
horizon is located at r = r + = r g (l + ^/l — (a/M ) 2 )). For 
a black hole with angular momentum J = aGM/c, a/M is 
the dimensionless Kerr pa rameter with — 1 ^ a/M ^ 1. 

The notation follows iMisner et al.l l|l973h and the sig- 
nature of the metric is (- ++. Tensor components are 

given in a coordinate basis. The components of the ten- 
sors of interest are given by j p for the metric, F^ v for 
the Faraday tensor, F for the dual of the Faraday, and 
T ,u/ for the stress-energy tensor. The determinant of the 
metric is given by \f—g = Det(gr M „). The field angular 
frequency is Qf = F tr /F r ^ = Fte/Fe^. The magnetic 

* it 

field can be written as B — F . The poloidal magne- 
tospheric structure is defined by the ^-component of the 
vector potential (A$). The current system is defined by 
the current density (J) and the polar enclosed poloidal 
current (B^, = F^ t ). The electr omagnetic lumi n osity is 
L = -2tt r p der (EM) t r 2 sing . See iGammie et all (l2003all ; 
McKinnev fc Gammiel l|2004h ; IMcKinnevI l|2004 l2005tl E 
2006al ) for details on this standard notation. Kerr-Schild 
coordinates are used t o avoid spurious reflections off the 
inner-radial boundary IjMcKinnev fc Gammie 2002). 



2 GRFFE MODEL WITH DISK 
MAGNETOSPHERE 

We study time-dependent GRFFE numerical models under 
the assumption that the turbulent accretion disk is well- 
modelled by an equatorial current sheet that is treated 
as a rotat i ng con ductor. The GRFFE code described in 
IMcKinnevI (|2006ah is used to evolve the axisymmetric force- 
free equations of motion in the Kerr metric in Kerr-Schild 
coordinates. This code has als o been successfull y used to 
study pulsar magnetosp heres (IMcKinnevI 12006^ 1 and the 
split-monopole problem (|McKinnevl l2006ah . Th e computa- 
tional g rid geometry is chosen to be the same as in McKinnev 
(2006c) with an inner radius inside the horizon and an outer 
radius at 10 3 r s . This grid focuses the resolution toward the 
polar axes at large radii in order to resolve the collimating 
jet. 

The resolution for all models is chosen to be 256 x 256, 
and the solutions are well-converged compared to low res- 
olution models except very close to the poles due to the 
coordinate singularity in spherical polar coordinates. 



2.1 Obtaining a Stationary Solution in GRFFE 

Steady state GRFFE numerical models with no discontinu- 
ities or surface currents above the disk surface are found 
by choosing boundary conditions determin ed by an analysis 
of the Grad- Shafranov equation (see, e.g.. lBogovalovlll997l ; 
lBeskirJll997h . For solutions that pass through a light cylin- 
der at some radius, one is required to fix the magnetic field 
component B e perpendicular to the disk and specify 2 other 
constraints at the disk, viz., = Ft^ = 0, and Qf = fidisk, 
the disk angular velocity. Note that one cannot leave both 
B e and B T unconstrained, as this allows arbitrary redistri- 
bution of the currents in the disk. Doing so gives too much 
freedom and the solution reverts to the split-monopole so- 
lution. No boundary condition is applied on the black hole. 

For axisymmetric, stationary solutions, the frozen-in 
condition of ideal MHD implies that the field line velocity 
if is completely determined by the field B l and field rota - 
tion frequency Qf (see equation 46 in IMcKinnevI l2006al ). 
Thus, during the simulation the 3-velocity in the equatorial 
plane is set to agree with this condition. Such a 3-velocity 
is generally time-like for points between the inner and outer 
light "cylinders," but outside this region the 3-velocity can 
be space-like and so unphysical. In the event that the 3- 
velocity is space-like, the Lorentz factor is const rained to a 
fixed la rge value as described in section 2.5 of IMcKinnevI 
(1200681 ). This safety feature is necessary to handle the vi- 
olent evolution seen early in the simulations, but it is not 
activated after the solution approaches a stationary state. 

In the work described here, the initial conditions of 
the simulations are chosen to correspond to a non-rotating, 
current-free (except at the equator) force-free solution with 
B^ = 0. For example, for the paraboloidal models, BZ's 
paraboloidal solution with a/M = is used to set the ini- 
tial cond itions. For other power-law c urrent profiles (e.g., 
v — 3/4) , Mc Kinnev fc Naravanl (|2006l ) describe how to cal- 
culate analytic solutions to the non-rotating problem. Dur- 
ing the time evolution, B 6 at the disk is held fixed at its ini- 
tial value, and the field angular velocity Qf is set to a phys- 
ical ly motivated profile as sugge sted by the GRMHD mod- 
els (|McKinnev fc Naravanll2006l '). For models with rotating 
disks and/or rotating black holes, the initial non-rotating 
state is far from the steady state solution, so the model un- 
dergoes violent non-stationary evolution. Eventually, how- 
ever, it relaxes to a steady-state solution. All solutions thus 
found are necessarily stable to Eulerian axisymmetric per- 
turbations. 



2.2 Modelling the Transition Region 

The transition region between the disk and the black hole 
must be treated carefully to avoid undesirable discontinu- 
ities in the magnetosphere away from the equatorial disk. 
We now discuss how the vector potential (A^, which is used 
to obtain B ) and the field angular frequency (Qf) are cho- 
sen. 



2.2.1 Softened Vector Potential 

If one uses an arbitrary vector potential to set the profile 
of B (r) at the disk, then this boundary condition may not 
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always allow for a stationary solution or one without dis- 
continuities. The problem is that we calculate B e from the 
vector potential (A,/,) for an idealized problem with a/M = 
or even M = 0, since these solutions are easy to obtain an- 
alytically. The solutions typically consist of power-law pro- 
files. However, for general M/0 and a/M ^ 0, the radial 
eigenfunctions must change character near the black hole in 
such a way that the solution for the vector potential must 
become nearly constant close to the horizon. 

This property of the force-free equations near the 
horizon is related to the fact that solutions have a min- 
imal energy if the black hole has a minimal t angen- 
tial field on the horizon iFrhorne fc MacDonaldl Il982al ; 
iMacDonald fc Thornd fl982bl T For example, co mpared to 



the p ure self-similar paraboloidal solution of iBlandfordl 
(1 19761 ). BZ's paraboloidal solution for a disk around a black 
hole effectively softens the vector potential near the black 
hole to a monopolar form in order to account for the pres- 
ence of the black hole. The characteristic radius where the 
change occurs in the case of a non-rotating hole is r ~ 2r g . 
For rapidly rotating black holes the horizon becomes smaller 
but the ergospheric radius in the equatorial plane remains 
at r = 2r 9 . So we expect generically that the softening of 
the potential will occur inside r ~ 2r g for all spins. 

In our work, we take the analytic vector potential so- 
lution corresponding to M = 0, a/M = and soften it by 
replacing the radial coordinate r as follows, 



A 2 , 2 

i r + T* trans j 



1/2 



(2) 



where r tra nsB is the transition radius. This allows the solu- 
tion to become monopolar near the horizon, thus modeling 
the effects of the mass and spin of the black hole. Usually we 
set rtransB = 2r+ , but none of the results depend sensitively 
on this choice, as we have confirmed by trying various values 
in the range r tra nsB = (1.5 — A)r g . 



2.2.2 Modelling the Angular Velocity 

The simplest force-free solutions have no discontinuities 
except as defined by the boundary conditions. In time- 
dependent force-free simulations, discontinuities can appear 
near the boundary if no simple force-free solution exists 
with the chosen boundary condition. This is something to 
be avoided. Therefore, in order to guarantee a simple force- 
free solution, the boundary conditions must sometimes be 
modified to match smoothly onto the magnetosphere. Even 
in quasi-analytic work it is useful to avoid discontinuities in 
the transi tion from the bounda ry conditions to the magne- 
tosphere i|Uzdensky||2003 . 120051 ) . 

In our models, the transition region where the disk 
meets the horizon requires careful treatment of the bound- 
ary condition in the disk. This is where the field line rota- 
tion switches from the profile set by the disk to the value 
required by the black hole space-time. We have tried two 
different models of the transition region. In one model, we 
take the field angular velocity profile to be given by 



Q. F = 



r <r t 



(l - { r -±yj + fi D BH (I^) nt , r > r transn . 

(3) 

The quantity Odbh denotes the value of Qf on the black 



hole horizon. We choose n t so that there is a smooth transi- 
tion. This model only applies for a/M ^ 0, and we find that 
values of r tr an S n in the range from r+ to 2r+ give acceptable 
results. 

A second transition model is chosen to enforce Qf to 
be strictly constant inside the horizon and strictly fixed to 
the disk profile beyond some transition radius. This model 
is given by 



{^dbh, r < r<), 

fioBH + A(r - r ) + B(r - r ) 3 , r ^ r < r trans n, 
f^disk, r > r t ransH 1 

(4) 

where J7dbh again denotes the value of SIf at the black hole 
horizon. The values of A and B are chosen such that Q,f and 
dQp/dr are continuous at ro and r tra nsn, where both radii 
are arbitrarily chosen. This model is suitable for any value 
of a/M. 

As in BZ77, we assume that on the horizon (r = r+) the 
angular frequency of field lines that thread the disk match 
onto some fraction of the black hole angular frequency 



0/ 



a 

2r~I 



(•») 



This matching of Qf between the disk and the magneto- 
sphere of the black hole enforces the condition that no extra 
surface currents (i.e., jumps in the field) are present that 
would have been created by a jump in the field angular 
frequency at the interface between the disk and black hole 
magnetosphere. To find force-free solutions with no disconti- 
nuities, we iteratively repeat the simulations with improved 
values of Odbh until the solution has a minimal discontinu- 
ity. 



3 BZ PARABOLOIDAL SOLUTION 

In this section the GRFFE equations of motion are self- 
consistently evolved to find the stationary force-free magne- 
tospheric solution corresponding to the toroidal current dis- 
tribut ion of the paraboloidal solution of Blandfor d fc Znaiekl 
|l977l ). viz., 



dr 



C 



(6) 



with v = 1. 

In the absence of rotation, the vector potential for the 
paraboloidal solution takes the form 



(para) 



A 



+g(r,0) 

+g(r,n-8) 



< vr/2 
>tt/2, 



C 



[r/_+2M/ + (l-ln/+ 



(7) 



(8) 



where /+ = 1 + cos 9, /_ = 1 — cos 8, and M is the mass of 
the black hole. As explained earlier, this non-rotating so- 
lution is used as the initial conditions for rotating mod- 
els, i.e., we obtain the boundary condition on the field at 
the equatorial plane from the non-rotating solution, viz., 
B e = —Afi,r/y?—g, and we hold this fixed for all time. We 
set C = 1 throughout the paper, and this normalizes all 
energy densities and field strengths. 
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Figure 1. Field angular frequency per unit black hole angular 
frequency (Qf/^h) vs. 6 in radians on the horizon (r = r g , solid 
line) and at two other radii (r = HSCO ~ 5.997r 9 and r = 10r 9 , 
dotted lines). The results correspond to a GRFFE simulation of 
the paraboloidal solution (y = 1) at time t = 1200t 9 for a non- 
rotating disk and a slowly spinning black hole with a/M = 0.001. 
Also shown are the Blandford-Znajek monopole solution (long- 
dashed line) and paraboloidal solution (short-dashed line) on the 
horizon. Note the excellent agreement between the GRFFE nu- 
merical solution and BZ's analytical model. 

The BZ paraboloidal solution is b ased on the non - 
gravitational solution found initially by iBlandfordl (1 19761 ). 
BZ hypothesized that the solution for the poloidal field is 
independent of the velocit y profi le of the disk (see eq. 3.10 
in that paper). Blandford] (|l976l ) noted in section 3iii that 
gravity will add corrections but, in fact, even without grav- 
ity there are corrections for any significant amount of disk 
rotation. Keplerian disks in particular have substantially dif- 
ferent solutions from the non-rotating case since the Keple- 
rian speed is a non-negligible fraction of the speed of light 
near the black hole . In th e discussion that follows, we show 
that the IBlandfordl (| 19761 ) and BZ solution for Qf and the 
energy output only applies when the disk is slowly rotating. 
We then calculate the correct solution for Qf and the power 
output for a Keplerian disk for the case of both slowly and 
rapidly rotating black holes. 

We model the transition region using equation @ with 
r = r+ and rtfansn = 2r+. We have also tried r tra n S n = 
1.2r+ and we find only small quantitative differences, though 
there are additional surface currents in the equatorial re- 
gions near the horizon because of the sharper transition in 
the Qf profile. 



3.1 Slowly Rotating BH and Non-Rotating Disk 

The first numerical model we describe has a/M — 0.001, 
i.e., a very slowly rotating black hole, and we assume a non- 
rotating disk, Sldisk = 0. We choose Odbh = 0. 27f2.tr as in 
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Figure 2. Angular density of the electromagnetic power 
output (dP/d.9) vs. 9 in radians at four radii, r = 
{?■-)-, TiscOi 10 r ff; 10 3 r 9 }. The results are for the same GRFFE 
simulation shown in Figure [T] The dashed line shows the analytic 
solution for the power output at the horizon according to the BZ 
paraboloidal model. Note the excellent agreement between the 
numerical GRFFE result and BZ's analytical model. 

the BZ solution. Figure [1] shows Qf/Qh for the converged 
steady state solution at time t = 1200t 9 at three different 
radii: the horizon at r = r+, the inner-most stable circular 
orbit (ISCO) at r — nsco, and r = 10r 9 . For comparison, 
the monopole and paraboloidal solutions on the horizon are 
also shown. 

The profile of Qf/Qh from the numerical solution is 
very close to the analytical BZ solution. The small offset is 
due to our model of the transition region, where one would 
expect to obtain BZ's solution only if allowing a disconti- 
nuity at the horizon as in their solution. There is a small 
residual jump near the equator in Qf due to not choosing 
the most optimal value for Qf on the horizon. Also, close 
to the poles Qf/Qh deviates from the expected value of 
Qf/Qh ~ 1/2 due to the coordinate singularity, though 
this feature has a negligible impact on the solution at large 
radii due to the collimation of field lines toward the polar 
axis. 

Figure [2] shows the angular density of the electromag- 
netic power output 

where _r (i5M) [ is written in a coordinate basis in Boyer- 
Lindquist or Kerr-Schild coordinates. The total integrated 
electromagnetic power at radius r is obtained by integrating 
over angle, 

r* rlP 

P= ^f-smOde. (10) 
Jo do 

The integrated electromagnetic power output at each 
of the four radii, r — {r+, nsco, 10r g , 10 3 r 9 }, is P = 
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{1.4, 1.6, 1.6, 1.6} x 10~ 6 , respectively, in units in which the 
radial field on the horizon at the poles is unity; that is, 
we scale the power as (dP/d9)/ Bq and P/Bq, where Bo is 
the radial field strength on the horizon at the polar axis 
at the end of the simulation. Normalized in the same way, 
BZ's paraboloi dal solution gives a power o f P ~ 1.2 x 10~ 6 
on the horizon (Bla ndford fc Z naick 1977). Clearly, our nu- 
merical solution is in good agreement. The small differences 
arise because our numerical solution for flp is slightly larger 
and the power on the horizon scales as P oc Qf(£If — 
We find that along field lines the value of Qf is constant 
to within 5%, as required for an axisymmetric, stationary 
solution. 

Figure [2] shows that the numerical profile for the angu- 
lar power density agrees well with the analytical BZ solu- 
tion. The slight difference is due to the introduction of the 
transition region between the disk and black hole. Note that 
the power becomes progressively more collimated with in- 
creasing radius. The peak in the power density lies at only 
9 « 10° by a radius of r = 10 3 r 9 . 

One can also compute the avera ge efficiency (e) over 



0.6 



the ho rizon, following equation 7.( 
(|1977| ): 



Blandford Znaiekl 



f n 
Jo 



' dP\ 

\ de J 



(11) 



Effectively, this is an estimate of Qf/CIh, weighted by the 
power density. This quantity has a maximum value of 50% 
for the monopolar solution. Our numerical solution gives an 
efficiency of 39%, which is similar to the value of 38% for 
the analytical BZ paraboloidal solution. 

In summary, the GRFFE numerical solution for the 
force-free magnetosphere surrounding a slowly-spinning 
black hole and a non-rotating disk agrees well with the an- 
alytical solution of BZ. 




Figure 3. Field angular frequency per unit black hole angular 
frequency (Qp /Qh) vs - in radians on the horizon (r = r g , solid 
line) and at two other radii (r = ngco ~ 2.321r 9 and r = 10r 9 , 
dotted lines). The results correspond to the paraboloidal solu- 
tion (y = 1) at time t = 1200f g for a Keplerian disk around 
a rapidly spinning black hole with a/M = 0.9. Also shown are 
the Blandford-Znajek monopole solution (long-dashed line) and 
paraboloidal solution (short-dashed line) on the horizon. The 
agreement between the numerical GRFFE model and the ana- 
lytical BZ model is not as good as for a slowly rotating black hole 
(Figure [T}, but the deviations are still fairly mild. 



of the transition region, we find additional surface currents 
in the magnetosphere above the disk. 



3.2 Slowly Rotating Black Hole and Keplerian 
Disk 

We now describe a second numerical GRFFE model, which 
is identical to the previous model except now the disk rotates 
at the Keplerian frequency, fidisk = fife = 1/ (r 3//2 + a) in 
Boyer-Lindquist coordinates. Because the Keplerian speed 
near the black hole is much larger than the black hole 
rotation speed, it leads to significant changes in the field 
around the horizon. The transition from the Keplerian disk 
to the relatively slow black hole speed forces the elec- 
tromagnetic energy to be accreted rather than extracted 
from the black hole. The integrated electromagnetic power 
output at the radii r — {r + , nsco, Wr g , 10 3 r 9 } is P = 
{—0.001,0.36, 1.42,2.04}, respectively, in the same units as 
before. For such slowly spinning black holes with rapidly 
rotating disks, the power output on the horizon depends 
sensitively on how we model the transition region between 
the disk and black hole. In particular, the negative power 
on the horizon results from a complicated set of currents in 
the magnetosphere that provide a black hole-disk connec- 
tion. The substantial power at large radii primarily reflects 
the power output of the rotating disk and there is very little 
power from the black hole itself. For all attempted models 



3.3 Rapidly Rotating BH and Keplerian Disk 

Finally, we consider a rapidly rotating black hole with 
a/M = 0.9, surrounded by a Keplerian disk (f2di s k = Ok). 
If we directly use BZ's paraboloidal solution for a slowly 
rotating BH to setup the disk boundary condition for the 
currents, only for a/M <C 1 were we able to find a solution 
without discontinuities near the horizon. In order to account 
for the additional effects of rapid spin, we had to soften the 
vector potential by replacing the radial coordinate as de- 
scribed in equation ([2]). This was not required for slowly 
spinning models. 

With this model of the transition region we were able 
to find a solution with Odbh = 0. 32f2.tr, which is a slightly 
larger value than for BZ's paraboloidal solution for a slowly 
spinning black hole. Figure [3] shows Qf/Qh at t = 1200i 9 
at three different radii and also shows the monopole and 
paraboloidal solutions on the horizon for comparison. 

Figure|4]shows the angular density of the power output. 
The peak in the power density occurs at a half-opening angle 
of 9°-14° at a radius of r = 10 3 r 9 . The integrated electro- 
magnetic power output at radii r — {r+,risco, 10r 9 , 10 3 r 9 } 
is P = {0.48,0.57, 1.39, 1.77}, respectively. Clearly, there is 
a non-negligible power output from the black hole and the 
transition region between the horizon and the ISCO. There 



© 2006 RAS, MNRAS 000.1111201 



Disk-Jet Coupling in Black Hole Accretion Systems II: FFE Models 7 




Figure 4. Angular density of the electromagnetic power 
output (dP/d8) vs. in radians at four radii, r = 
{f+,nscO) 10 r s; 10 3 r g }, for the GRFFE solution shown in Fig- 
ure \3\ The power is collimated within a half-opening angle of 
9°-14° by r = 10 3 r s . 



is substantial power from the disk as well. From equation 
the average efficiency is 41% for this model, which has 
a rapidly rotating black hole. It is slightly higher than in the 
original BZ model for a slowly spinning black hole. 

Figure [5] shows the structure of the magnetosphere 
once the solution has reached a stationary state. As in the 
paraboloidal solution of Blandford-Znajek, the field geome- 
try becomes monopolar near the horizon. This general rela- 
tivistic effect becomes prominent once the field lines thread 
the ergosphere. 

Now we consider the radial dependence of certain quan- 
tities along a field line starting at 9 w 60° on the horizon 
out to large radii. The half-opening angle of the field line is 
found to follow 



54° 



2.8rv 



(12) 



which is close to the scaling 6 j oc r" 1 ^ 2 for the non-rotatin g 
v = 1 paraboloidal solution (Mc Kinnev fc Naravanl feOQG). 
The field has apparently slightly decollimated compared to 
the non-rotating case, which is also clear from Figure [6] 
where the field lines in the rotating model are seen to be 
slightly decollimated with respect to the nonrotating solu- 
tion. Thus it appears that the decollimating centrifugal force 
associated with rotation is somewhat stronger than the col- 
limating hoop-stresses associated with the toroidal field gen- 
erated by rotation. 

For a force-free solution, the minimum Lorentz factor 
with which observer at infinity would see particles moving 
is given by the drift speed. At large radii this minimum 
Lorentz factor is well-fitted by 



r = u y/-g t t ~ 1.5 



42.5r c 



(13) 



Figure 5. For the paraboloidal (y = 1) model with a/M = 0.9 
and a Keplerian disk, this shows the poloidal field geometry given 
by contours of the vector potential (A^). Also shown, from in- 
ner to outer radius, are the inner-radial computational boundary, 
horizon at r = 1 + vl — a 2 , the inner "light cylinder" corre- 
sponding to the Alfven surface for ingoing Alfven waves, and the 
ergosphere at r = 1 + y/1 — (a cos 0) 2 . Notice that the field ge- 
ometry becomes monopolar near the horizon as in the original 
Blandford-Znajek solution. 



which is similar to the scaling Y oc r 1 ' 2 found by 
iBeskin fc Nokhrinal l|2006l ) for paraboloidal MHD models. 
This suggests that the force-free flow may be a good model 
of the acceleration regime in GRMHD models, at least at 
those radii where the magnetic energy is much larger than 
the kinetic energy. 

The orthonormal toroidal field at large radii follows 



B^ w B^^fg^ w 0.26 I — 



and the pitch angle at large radii follows 



"pitch = tan 



tan 



(14) 



(15) 



where the speed of light (c) has been temporarily reintro- 
duced. This is consistent with an orthonormal radial field of 
the form 



B r w -- 



-B 4 



(16) 



which is a generic feature of force-free solutions with rotat- 
ing field lines (see equation ()A8|I in Appendix|X]). Studies of 
stability suggest that force-free fields that follow the above 
behavior are stable to th e nonaxisymmetric kink instability 
l|Tomimatsu et all 1200 ll ). though our axisymmetric simula- 
tions cannot test this. 
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R cVCM 

Figure 6. Field lines (constants of As) for the a/M = 0.9 
paraboloidal [u = 1) model at t = (initial state, non-rotating 
solution, dotted lines) and t = 1.2 X 10 3 t 9 (final time, converged 
rotating solution, solid lines). The field is seen to decollimatc 
slightly in the presence of rotation, but the effect is quite mild. 
Slower and more rapidly rotating black holes show similar results. 

3.4 Summary of Paraboloidal Solutions 

The original BZ paraboloidal solution could not describe the 
transition region between the accretion disk and black hole, 
it does not account for modification of the field direction for 
disks with Keplerian rotation profiles, and it only applies to 
slowly spinning black holes. Our numerical work has elimi- 
nated these restrictions. For rapidly rotating black holes, we 
find an enhanced averaged efficiency for the power from the 
black hole. We find that power output from the transition 
region between the disk and black hole contributes a signif- 
icant fraction of the total power output from the system. 
We find that for slowly spinning black holes the transition 
region is important in determining force balance and the 
power output on the horizon. 



4 GRMHD-MOTIVATED v = 3/4 SOLUTION 

Since GRMHD numerical models pro duce disk currents con- 
sisten t with a v = 3/4 distribution l|McKinnev fc Naravanl 
2006), we now focus on force- free solutions with this scal- 
ing for the current. In this section we describe the model of 
the toroidal current and the field rotation frequency as de- 
duced from GRMHD simulations, and we study correspond- 
ing GRFFE numerical models to determine the power and 
collimation of the Poynting-dominated jet and electromag- 
netic disk wind. 

4.1 Current Model 

The initial conditions are chosen to be the non-rotating 
flat-space (M = a/M = 0) f = 3/4 model described in 



iMcKinnev fc Naravanl l|200 6:) . The vector potential is given 
by 

A^ = Cr " | sin 6» | Pj_! (cos 6»), (17) 

where P™ is the associated Legendre function of the first 
kind. For the case of interest, viz., v — 3/4, the angular 
solution takes the form P" 1 with I = —1/4 and m — 1. 
We set C — 1 throughout the paper, and this normalizes 
all energy densities and field strengths. As in the case of 
the paraboloidal models, this vector potential is used to ob- 
tain the boundary condition on B at the disk, which is 
then held fixed for all time. We modify the self-similar solu- 
tion to account for the transition region near the black hole 
via the mapping given in equation ((2]). The main effect of 
this smoothing operation near the black hole is the intro- 
d uction of a split-monopole fie ld ne ar the center. As f ound 
in IMcKinnev fc Gammiel (|2004l ) and lMcKinnevI (l2006d ). the 
field is indeed nearly monopolar across the horizon in the 
Poynting-dominated jet region. This is consistent with the 
fact that all force-free solutions with M/0 tend to become 
monopolar near the horizon. 

4.2 Field Rotation Model 

We use GRFFE numerical models to investigate the two 
models of SIf described in equations © and @. For the 
latter, we choose f2disk = ro = 2r + and r tra nsn = 3r + , 
though we have found that choosing ro = r+ and rtranso = 
2r+ leads to only small quantitative differences. To fit the 
ti me-averaged GRMHD numeri cal model shown in Figure 8 
in IMcKinnev fc Naravanl (|2006l ). one could choose ^dbh = 
0.6Qh- For general a/M, a good fit is provided by Hdbh ~ 
S7h/2 for a /M > 0.4 and Odbh = for slower black 
hole spins (|McKinnev fc Naravanl 120061) . Note that Q F is 
the same in Kerr-Schild and Boyer-Lindquist coordinates, 
and Q is also the same at the equator. Thus, the above 
prescription for f2f? in the equatorial disk is valid for both 
coordinate systems. 

4.3 Rotating Black Hole and a v — 3/4 Keplerian 
Disk 

We consider a GRFFE model with a black hole spin of 
a/M = 0.9 and an equatorial current sheet with & v — 3/4 
power law current distribution. Our fiducial model for the 
transition region, called "TModel 1," smooths the field vec- 
tor potential using equation © with ftransB = 2r+, and 
models the field rotation via equation ((3} with Odbh /fijj = 
0.32 and r tra nsn = 2r+. 

For this transition model we seek to find force-free solu- 
tions with no discontinuities. If one chooses Odbh = 0. 6f2.tr 
as in the GRMHD simulations, discontinuities appear above 
the equatorial plane. This is expected since the disk is not 
force-free. The value Odbh = 0. 32f2.fr was chosen to give a 
smooth transition from the disk to the black hole so that 
there are negligible discontinuities near the accretion disk - 
black hole interface that would correspond to extra toroidal 
currents in the magnetosphere. Other transition models are 
considered in section 14.51 such as one where we do choose 
Odbh = 0.6f2j/ and allow discontinuities to be present in 
the magnetosphere. 
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Figure 7. Field lines (black lines), enclosed poloidal current (Bj,, 
where B& yjg$$, red positive, blue negative), and light 

"cylinder" (thick black line). The results correspond to a GRFFE 
model with v = 3/4 at time t = 1200t ff for a rotating disk around 
a spinning black hole with a/M = 0.9. The panel on the left 
shows the full computational domain and that on the right shows 
a closeup of the central region. The black hole and the disk to- 
gether generate a Poynting outflow in the form of a jet surrounded 
by a wind, with power focused along the polar axis. Small artifacts 
at the outer radial edge in the left panel are the result of inter- 
polation and reflections of the Poynting jet off the outer radial 
boundary. 

Starting with the non-rotating force- free solution as ini- 
tial conditions, the above model was simulated with the 
GRFFE code until it reached steady state (t = 1200i 9 ). Fig- 
ure [7] shows the poloidal structure of field lines (also lines 
of constant A$), the enclosed poloidal current from the pole 
(Btf,), and the light "cylinders" (only the outer light cylin- 
der is clearly visible) in the final state. The poloidal field 
geometry has changed little from the initial solution, show- 
ing that rotation has a negligible effect. The polar enclosed 
poloidal current (B^) is shown in color, where the radial 
energy flux is given by — QfB^B 7 ^ . Note that for sta- 
tionary, axisymmetric flows i?0 and SIf are constant along 
field lines. 

The field lines in the v = 3/4 solution show less poloidal 
collimation than the paraboloidal solution but more than the 
monopole solution (which has no collimation at all) . As com- 
pared to the GRMHD solution, the pressure support of the 
corona has been replaced by the stable force-free field lines, 
but otherwise the funnel region is quantitatively similar to 
the GRMHD solution. 

Figure [8] shows the structure of the magnetosphere 
once the solution has reached a stationary state. As in the 
paraboloidal solution described previously, the field geome- 
try becomes slightly monopolar near the horizon. This gen- 
eral relativistic effect becomes important once the field lines 
thread the ergosphere. 

Figure[9]diagnoses the poloidal current flow and angular 




R c 2 /GM 



Figure 8. Similar to Figure [5] but for the GRMHD-motivated, 
a/M = 0.9, Keplerian disk model with a v = 3/4 current distribu- 
tion. Notice that the field geometry becomes slightly monopolar 
near the horizon. 

rotation of field lines in the force-free magnetosphere at a 
radius of r = 20. The value of B^ is the amount of poloidal 
current enclosed away from the polar axis down to the given 
value of A$, where A$ = on the polar axis. This plot shows 
that there is an increase in current away from the polar axis 
due to the jet coming from the black hole, and then there 
is a drop in current as the disk supplies the return current. 
At large radii (large A$) the current vanishes. The current 
closes at r = and r = oo in a non-singular manner, unlike 
in the paraboloidal case where B^, is finite at large radii 
so that the radially integrated poloidal current is infinite. 
The small-scale oscillation in B^ is located at the black hole 
jet-disk interface, and the detailed behavior of this region 
depends on the black hole spin and disk rotation profile. 
This current structure is c omparable to that of neutron star 
magnetospheres (see, e.g. McKinncy 2006b). The plot also 
shows Qf /SIh as a function of A$. This shows how the black 
hole solution at small A$ matches onto the Keplerian disk 
at large A$. 

4-3.1 Angular Dependence within Jet 

Figure [TU] shows Qf/Qh for the v = 3/4 solution at t = 
1200t 9 at three different radii and also shows the monopole 
and paraboloidal solutions on the horizon. As expected, 
the v = 3/4 solution has a field rotation at the horizon 
(£If/£Ih ~ 0.32) that is close to the paraboloidal solution 
(v = 1, Qf/^Ih ~ 0.32 for a/M = 0.9) and lower than 
the monopole solution (u — 0, Hf/Hh ~ 0.5). The profile 
of SIf/Qh is similar to the GRMHD mode ls show n in Fig- 
ure 8 (left panel) in iMcKinnev fc Gammiei (|2004| ). In their 
plot the value of £If/Qh ~ 0.32 near the transition between 
the force-free region and the accretion disk region. This is 
further evidence that the v = 3/4 GRFFE model gives a 
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Figure 9. Enclosed poloidal current (B$, where B$ « B^^Tg^, 
solid line with crosses) and field angular frequency per unit black 
hole angular frequency (Qp/Qx, dashed line) vs. the vector po- 
tential [Aj, = \P, also called the flux function) at spherical polar 
radius r = 20. This shows that the poloidal current increases in 
the jet and then decreases at larger radii in the disk that con- 
tributes to a line return current. There is no enclosed current at 
r = oo unlike in the paraboloidal case. This also shows how the 
rotation frequency of the field lines behaves near the black hole 
and how it connects onto the Keplerian disk at large radii (large 



Figure 10. Field angular frequency per unit black hole angular 
frequency (Qp /Qh) vs - 6 m radians on the horizon (r = r g , solid 
line) and at two other radii (r = rigccj ~ 2.3209r 9 and r = 10r 9 , 
dotted lines). The results correspond to a GRFFE model with 
v = 3/4 at time t = 1200t 9 for a rotating disk around a spinning 
black hole with a/M = 0.9. Also shown are the Blandford-Znajek 
monopole solution (long-dashed line) and paraboloidal solution 
(short-dashed line) on the horizon. Note that the field angular 
frequency at the horizon is roughly similar for the paraboloidal 
(v = 1) and v = 3/4 models. 



solution that is consistent with the Poynting-dominated jet 
found in full GRMHD. Close to the poles flp/Qn deviates 
from the expected value of Qp/Qh ~ 1/2 due to the coordi- 
nate singularity. However, this feature has negligible impact 
on the solution at large radii due to the collimation of field 
lines toward the polar axis. 

Figure [11] shows the angular density of the power out- 
put. The v = 3/4 disk provides most of its energy inside 
r = 100r 9 . The integrated electromagnetic power output at 
r = {r+,nsco,10r 9 ,10 3 r 9 } is P = {0.47,0.62,1.50,1.76}, 
respectively. The figure shows that the peak in the power 
is collimated to within a half-opening angle of 15°-26° by 
r — 10 3 r 9 . Using equation (|11[) . we calculate an average 
efficiency of 41% for this model, which is the same as the 
equivalent paraboloidal model. 

One can compare the electromagnetic power outputs in 
the GRFFE model to those found in the equivale nt GRMHD 
model discussed in lMcKinnev fc Naraya n (2006). The pow- 
ers are comparable at the horizon. However, in the GRMHD 
model, the electromagnetic power from the disk is lost to the 
matter and only the power from the black hole survives at 
large radii, whereas in the GRFFE model all the electromag- 
netic power output of the black hole and the disk reaches 
large radii. 
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Figure 11. Angular density of the electromagnetic power 
output (dP/d9) vs. 8 in radians at four radii, r = 
{r+,rjgco! 10r 9 , 10 3 r 9 }, for the GRFFE solution shown in Fig- 
ure [10] The power is seen to be collimated within a half-opening 
angle of 15°-26° by r = 10 3 r 9 . The collimation is significantly 
less than in the otherwise similar paraboloidal (y = 1) model. 
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Figure 12. Variation with radius of various quantities in the converged GRFFE solution with v = 3/4 and a rotating disk around a 
spinning black hole with a/M = 0.9. The solid lines in the panels show the half-opening angle of the jet in degrees (0j), the Lorentz factor 
(r), the comoving energy density (fe 2 /2), the orthonormal radial velocity (u r ), the orthonormal toroidal field strength (-B^), and the 
pitch angle in radians (ctpitch)- The features near the outer radial boundary are due to reflections of initial transients off the boundary. 
Except for the panel with a p itchi the dashed lines in the other panels show power- law fits to the numerical results. The dotted line in the 
u r panel shows the orthono rmal angular veloc i ty . The dashed line in the a p it c h panel shows the limit below which the kink instability 
will be present according to lTom imatsu et al. (2001). It indicates that the flow is marginally stable at large radii. The other panels show 
that the solution has a power-law behavior at large radii. 



4-3.2 Radial Dependence within Jet 

Figure [121 shows the jet half-opening angle (0j), the Lorentz 
factor (r = i^yf—gtt), the comoving energy density (b 2 /2), 
the approximate orthonormal absolute value of radial and 
4> velocities (\u r \ = \u r \^Jg rr and = u^^g^), the or- 
thonormal toroidal field strength (B^ = B^^/g^), and the 



pitch angle within the jet (a p itch = tan _1 (B r /B^)). The fig- 
ure also shows power-law fits to some quantities. This figure 
can be compared to panels in figures 7 and 8 in iMcKinnevI 
l|2006ch . 

The opening angle of the jet for the field line starting 
at 0j ~ 57° has a power-law dependence at large radii that 
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follows 



0< « 57° 



2.8rv 



(18) 



comp ared to 9j oc r~ 0,375 for the non-rotating v = 3/4 solu- 
tion l|McKinnev fc Naravanll2006fl . This shows that the field 
has slightly decollimated compared to the non-rotating solu- 
tion, but the effect is weak. The power-law dependence for 9j 
is strikingly similar to that found for the opening angle of the 
je t in GRMHD simul ations given by equations (24) and (25) 
in iMcKinnevi (|2006cD . For example, the choice of Oj « 57° 
on the horizon for the field line corresponds to where the 
Po ynting-dom i nated jet starts in the GRMHD simulations 
of IMcKinnevi (2006a) an d corresponds to where the com- 
bined disk+corona scale height to radius ratio is H/R « 0.6. 



In this force-free mod el, 6j « 5° at r 



5 x 10 3 r 9 , a similar 



value as found in IMcKinnevi (|2006cT ) for the core of the jet 
at large radii. 

The GRMHD simulations and the present force-free 
model have an opening angle that is consistent up to the 
point where in the GRMHD model the jet becomes coni- 
cal due to dissipative processes. In the GRMHD simulations 
the power-law dependence slightly depends on the field line 
chosen, while in these GRFFE simulations the power-law 
dependence is the same for all field lines. 

The Lorentz factor at large radii is well-fitted by 



(19) 



42.57% 



which is consi stent with t he ma gnetic acceleration do- 
main found in IMcKinnevi l|2006ch . This behavior is also 
similar to the ana l ytic r esult of F oc r 1 / 2 obtained by 
Bcsk in fc Nokhrinal (|2006T ) for para boloidal MHD win ds. 

As in the GRMHD model of IMcKinnevi l|2006ch . the 
plasma angular velocity (u*) remains subrelativistic while 
the radial velocity (u r ) becomes relativistic. 

The comoving energy density follows 



2 



1.4 



(20) 



which is close to the scaling b 2 oc r -2,5 found for the non- 
rotating v — 3/4 solution. This implies that the black hole 
and disk rotation has only a small effect on the field energy. 
The orthonormal toroidal field follows 

B* « 0.29 (^) , (21) 

which is remarkab ly close to the GRMHD solution found in 
IMcKinnevi (|2006cT ) for the inner-radial region before dissipa- 
tion becomes important at r ~ 10 3 r 9 . 
The pitch angle is found to follow 



Opitch ~ tan 



(22) 



with the speed of light (c) temporarily reintroduced. As 
for the paraboloidal models, the orthonormal radial field at 
large radii is well-fitted by 

B* m — — — — B*. (23) 



Summarizing this section, there are two main results. 
First, the poloidal structure of the field in the rotating 



GRFFE models considered in this paper is rather similar 
to the fiel d structure found in th e non-r otating models dis- 
cussed in iMcKinnev fc Naravanl l|2006h . Since we showed 
previously that the latter models agree well with the poloidal 
field structure in the jet region of GRMHD models, the 
present models also agree well. Second, the structure of the 
toroidal field, along with various related quantities such as 
the pitch angle, Poynting flux, jet acceleration, Lorentz fac- 
tor, etc., obtained from the rotating GRFFE solution agree 
very well with the corresponding quantities in the GRMHD 
jet. Thus, the rotating GRFFE model appears to be an ex- 
cellent representation of all facets of jets. 



4.3.3 Stability 

The pitch angle of the magnetic field lines becomes pro- 
gressively smaller at large distances, becoming as small as 
Qpitch ~ 3° at r — 10 3 r 9 . This is consist e nt wit h the 
Poynting-dominated jet studied in IMcKinnevi (l2006ch . The 
small pitch-angle argues that the flow might be kink unsta- 
ble, since the Kruskal-Shafranov criterion for instability is 
abundantly satisfied. The criterion for instability is given by 



B <t> 

w 



~ V 



(24) 



where L « ^/gee is approximately the vertical length of the 
jet and R « ^fg^, is approximately the cylindrical radial 

exte nt of the jet. 

iTomimatsu et ah! l|200 ll ) considered the effect of rota- 
tion on stability and found that a Poynting-dominated jet is 
kink unstable only if both the Kruskal-Shafranov criterion 
and the additional condition 



R_ 

Rl 



(25) 



is also satisfied, where R l = c/Qf is the light cylinder radius 
|Tomimatsu et all 1200 ll ). Their result only strictly applies 
inside the light cylinder and for a uniform field, but it does 
suggest that rotation could stabilize a jet that might ap- 
pear unstable according to the Kruskal-Shafranov criterion. 
The lower right panel of Figure [T^] compares the quantity 
tan" 1 (Rl / \jg<j>$) with the pitch angle, and suggests that the 
flow is marginally stable to the kink instability at large radii. 
Note, however, that the kink instability is non-axisymmetric, 
whereas our simulations are axis ymmetric. Thus, 3D sim - 
ulations are needed to check the ITomimatsu et all (|200ll ) 
stability analysis. 

We can, however, study the stability of our numeri- 
cal force-free solutions to the axisymmetric pinch instability 
both inside and outside the light cylinder. For a/M = 0.9, 



line that threads the black hole (as given by equation (|18[) . 
the outer value of the cylindrical radius is R w 120r 9 . Thus 
such a field line is far beyond the light cylinder. Field lines 
that thread the outer disk have not yet passed through their 
light cylinder, so we focus on the stability of the solution 
close to the black hole. 

Once the solution became stationary by t ~ 10 3 t 9 , ran- 
dom perturbations with an amplitude of 10% were added. 
The Poynting-jet from the black hole and the electromag- 
netic wind from the disk were found to be stable to these 



© 2006 RAS, MNRAS OOO.ITH201 



Disk-Jet Coupling in Black Hole Accretion Systems II: FFE Models 13 



perturbations within and far beyond the light cylinder where 
the toroidal field dominated. Thus, we can state that the 
axisymmetric pinch mode does not grow, at least in an Eu- 
lerian sense (though it may still grow in a Lagrangian sense 
moving out with the fluid and may manifest itself at large 
radii). 

4.4 Dependence on Black Hole Spin 

Now we study the dependence of the power on the black 
hole spin for GRMHD-motivated GRFFE models. In the 
GRMHD simulations, the value of Qf/Qh in the disk de- 
pends strongly on whether the spin is larger or smaller than 
a/M ~ 0.4, thus a mapping to a force-free model is not 
straight-forward because such force-free models choose an 
approximately fixed value of Qf/Qh for the entire mag- 
netosphere above the equa torial plane. As described in 
IMcKinnev fc Naravanl < |2006h . GRMHD models with black 
hole spins a/M > 0.4 tend to have Qf/Qh ~ 1/2 on the 
horizon, and so using Qf/Qh = 0.32 in the GRFFE model- 
ing should lead to a solution close to the GRMHD solution. 
However, for a/M < 0.4, the GRMHD solution within the 
disk is found to follow Qf/Qh ~ 1 because the disk dom- 
inates the field all the way through the horizon. Since the 
power from the black hole is P oc Qf(Qf — QH)(B r ) 2 , then 
P ~ for a/M < 0.4. 

We already know that the magnetosphere is free from 
unwanted discontinuities only if one chooses Qf/Qh ~ 0.32 
on the horizon at the equator. One problem with mod- 
elling a/M < 0.4 black holes is that choosing the value 
of Qf/Qh ~ 1 would lead to discontinuities in the mag- 
netosphere. Even if one allowed discontinuities, the more 
important problem is that the field lines passing through 
the horizon above the equator would still settle on having 
Qf/Qh = 0.32. Hence, while the power through the disk 
would vanish exactly at the equator, it would not vanish just 
off the equator. A comparable GRMHD simulation with a 
thick disk would have P ~ over the entire part of the hori- 
zon occupied by the infalling disk material. In such models, 
as the disk thickness increases, the power (per unit (B r ) 2 ) 
drops because the disk takes over the magnetosphere. Thus, 
the disk thickness plays an important role in setting the total 
power output for a/M < 0.4. 

Despite these complications, in the following we employ 
the same model of Qf/Qh in the disk for all black hole spins, 
but we primarily focus on the models with a/M > 0.5 where 
the results are most representative of the GRMHD models 
for any disk thickness. We use equation (J5J to model the 
disk rotation and equation ([2]) to soften the vector potential 
in the transition region. 

Table Q] shows the electromagnetic power output, 
Lorentz factor, and half-opening angle for different black 
hole spins and at different radii. The power could be arbi- 
trarily normalized, but we use two specific normalizations 
that are convenient for comparisons with the results from 
GRMHD models. In the upper section of the table, the power 
has been normalized by (B r ) 2 on the horizon at the poles 
of the initial a/M = M = model with v = 3/4. Since the 
field strength at large radius is not significantly changed by 
the slow disk rotation or the rotation of the black hole, this 
is equivalent to a normalization by the disk field strength at 
large radii. If the accretion rate and field strength at large 



radii can be established for a GRMHD model, then this is a 
good normalization to use for comparison with the GRFFE 
model. The lower section of the table has the power normal- 
ized by (B r ) 2 on the horizon at the poles of the final station- 
ary model. This removes the changes in the field strength 
near the black hole due to the spin of the black hole. If the 
accretion rate at small radii is in force-balance with the field 
near the black hole, then this is a good normalization. How- 
ever, the field stre ngth near the hori zon is a strong function 
of black hole spin (|McKinnevll2005al ') and a mild function of 
the mass accretion rate, so this normalization should be con- 
sidered carefully when comparing results for different black 
hole spins. 

Regardless of the normalization, one interesting conclu- 
sion, from the table values for -P[r+] to P[r = 10 3 r 9 ], is that 
even at high black hole spins the disk provides no less than 
« 1/4 of the total power output. 

The Lorentz factor and 9j are independent of the nor- 
malization of course. The power shifts to small angles due 
to the emergence of the black hole component above a/M ~ 
0.5, although this black hole component is always present. 
For the a/M = 0.9 model, the field line that connects to 
the black hole at the equator eventually reaches 0j w 10° at 
r = 10 3 r 9 . The peak in the angular power density is set by 
the black hole power output for a/M > 0.9 and is otherwise 
set by the disk power output. 

If a disk were present with a sharply defined height 
to radius ratio of H/R w 0.3, as in the GRMHD 
numerical models, and the matter did not modify the 
poloidal field direction (whi ch is a good approximation; see 
IMcKinnev fc Naravanl 120061 ). then the field line would go 
from 9 « 72° on the horizon to 9 3 « 10° at r = 10 3 r 9 . In 
the presence of both an accretion disk and corona, as in the 
GRMHD numerical models, with a combined scale-height to 
radius ratio of H/R « 0.6, then 9 3 « 7° at r = 10 3 r 9 . Fi- 
nally, if we extend the calculation to r = 5 x 10 3 r 9 and again 
consider a disk+corona region blocking part of the power of 
the jet, then 9j w 5° there. This is the same jet collimation 
found for the co re of the jet in tim e-dependent GRMHD nu- 
merical models (|McKinnevll2006d\ In those simulations the 
jet had an extended low-energy wing out to 9j ~ 27° , which 
could be due to the effects of matter or due to a lack of 
transfield pressure support at large radii. Overall the open- 
ing angle of the jet and the power output measured in the 
GRMHD numerical models and the v = 3/4 force- free mod- 
els are in excellent agreement. A monopolar, paraboloidal, 
or cylindrical force-free model would not agree this well. 

Figures [T3l and [T4l show the initial and final configura- 
tions of the field lines for the a/M = 0.1 and a/M = 0.9375 
models, respectively. We see that the Keplerian rotation of 
the disk leads to some collimation of the disk field lines. On 
the other hand, the spin of the black hole leads to modest 
decollimation of the field lines near the poles. 



4.5 Dependence on the Transition Model 

Different models of the transition region were studied and 
Table [2] gives the same information as Table [T] for these dif- 
ferent models. 

The first model, called "TModel 1," has been already 
discussed in §§ I4.3l and l4.4l It uses equation ((3| for Qf with 
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Table 1. Electromagnetic Power, Lorentz factor, and half-opening angle of jet for different black hole spins and at a few radii for each 
model. The second to last column shows 8j at the location of the peak in the angular power, which includes a significant power from 
disk at low black hole spins. The last column shows 8j at a large radius for the field line that starts at 6j = 57° on the horizon, which 
begins the force-free region in GRMHD numerical models that have an accretion disk up to H/R s=s 0.3 and a corona up to H/R fa 0.6. 
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TModel 1 


2.097 


2.744 


6.662 


7.846 


15 


12°- 


-26° 


5° 


TModel 2 


1.672 


1.938 


3.983 


4.65 
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-26° 
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TModel 3 


10.64 


11.47 


20.82 


21.84 


15 


12°- 
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5° 



Normalized by Black Hole Field Strength 

TModel 1 0.4715 0.6168 1.498 1.764 15 12°-26° 5° 

TModel 2 0.4998 0.5794 1.191 1.39 15 15°-26° 5° 

TModel 3 0.5063 0.5456 0.9906 1.039 15 12°-26° 5° 



Table 2. Different models of the transition region are studied for a fixed black hole spin of a/M = 0.9. Otherwise similar to Table [T] 



^dbh/^jj = 0.32 and rtransfi = 2r + , and equation ((2)1 for 
the vector potential with rtransB = 2r+ . 



TModel 2 

The second model, called "TModel 2," uses equation Q for 
Q F with Q VB h/Qh = 0.32, r tra nsn = 3r + , and r = 2r + , 
and equation ([2]) with r tr ansB = 2r+. This model can be 
used to study any a/M and has a smooth transition from a 
Keplerian Hf at large radii to a fixed Qf near the black hole. 
TModel 2 leads to qualitatively and quantitatively similar 
results as TModel 1. 



TModel 3 

The final model, called "TModel 3," uses equation (|3]) for 
Qf with Qdbh/Qh = 0.59 and r tra nsn = 2r+, and again 
equation Q with rtransB = r+. This last model is chosen to 
match most closely to the GRMHD transition regardless of 
whether additional discontinuities are created in the mag- 
netosphere. 

Most features of this model are qualitatively similar, 
such as the field lines as shown in Figure Q3] and the power 
output distribution as shown in Figure 1111 There are two 
new qualitative results, however. First, there are additional 
current sheets in the magnetosphere near the disk inside the 
ergosphere. This is due to £If not matching between the 
disk and the black hole. These additional currents indicate 
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R cVCM 

Figure 13. Field linos for the a/M = 0.1 GRFFE model with 
v = 3/4 at t = (initial state, nonrotating solution, dotted lines) 
and t = 1.2 X 10 3 i 9 (final converged rotating solution, solid lines). 
The field lines threading the black hole show mild decollimation, 
as in the paraboloidal case, and the field lines from the outer 
regions of the disk show some collimation. 




R cVCM 

Figure 14. Similar to Figure IT31 but for a rapidly spinning black 
hole with a/M = 0.9375. There is again mild decollimation of the 
field lines threading the black hole, and some collimation of field 
lines from the disk. 



that no simple force-free solution can replace the matter- 
dominated disk. The second feature is that the field strength 
at the poles of the black hole is about 2.5 x more enhanced 
compared to the initial field strength. This field enhance- 
ment and the faster disk rotation leads to about 5x more 
power output from the black hole and about 2x more power 
from the disk if one normalizes the power by the disk field 
strength. This suggests that an accurate model of the tran- 
sition region, only possible in full GRMHD, is important in 
obtaining a quantitative estimate for the power output. 

Transition Model Summary 

In summary, the accreting matter plays a significant role in 
setting up the transition region, and no trivial normalization 
of the force-free results can provide an unbiased estimate 
of the power output. In GRMHD models, the disk+corona 
absorbs much of the power that reaches large distances in 
the GRFFE models. Force-free models with a thick wedge 
to replace the disk+corona may better reproduce the power 
lost to the disk and corona. 

However, we found that the qualitative behavior of the 
force-free Poynting-jet from the black hole does not depend 
on the details of the matter flow in the transition region be- 
tween the disk and black hole. For a good estimate of the 
power output that would correspond to the power output 
found in GRMHD numerical models, one can truncate the 
power output in the GRFFE model to include that power 
only above the disk+corona scale height in the funnel that 
contains the Poynting-dominated jet. For the GRMHD mod- 
els studied here, this angle away from the disk is H/R « 0.6. 
To convert the power output to physical units, one can use 
the dependence of the horizon magnetic field strength on the 
mass ac cretion rate and b lack hole spin found in GRMHD 
models l |McKinnej2005af ). 

Quantitative estimates of the Lorentz factor and open- 
ing angle are insensitive to the details of the transition re- 
gion. This is because with a v — 3/4 current distribution, 
the poloidal field direction does not significantly change for 
any changes in the disk or black hole angular frequency. 

4.6 Comparison of GRMHD and GRFFE 
a/M = 0.9375 Models 

Figure [15] shows the field geometry from a time-dependent 
GRMHD numerical model as shown in Figure 9 from 
iMcKinnev fc Naravanl ((2006) , but now overlayed with the 
a/M = 0.9375 v = 3/4 GRFFE stationary solution for the 
poloidal field. The funnel region, along which the Poynting- 
dominated jet in the GRMHD model emerges, has a poloidal 
field reasonably consistent with the v = 3/4 current sheet 
solution. 

Hoop-stresses can be discarded as important in help- 
ing to collimate the jet since the poloidal field geometry is 
decollvmated by the action of the black hole spin, leading to 
even larger differences between the GRMHD and GRFFE 
poloidal fields near the poles compar ed to the differences 
found in IMcKinnev fc Naravanl (|2006l ) where we compared 
the polar field from the GRMHD solution and the non- 
rotating v — 3/4 force-free solution. 

The field lines are more collimated in the GRMHD 
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Figure 15. Overlay of the time-averaged poloidal field from the 
GRMHD numerical model of an accretion disk with a rapidly 
spinning black hole with a/M = 0.9375 (black) and the poloidal 
field from the numerical GRFFE model with v = 3/4 and the 
same value of a/M (red). Only the portion of the v = 3/4 GRFFE 
solution that overlaps the funnel region of the GRMHD model is 
shown. Note the excellent agreement between the two models in 
the overlap region. However, the agreement is no better than with 
the non-rotating v = 3/4 solution shown in Figure 9 of McKinney 
Sz Narayan (2006). Thus, while the choice v = 3/4 is key to obtain 
a good match in the poloidal field structure, the rotation of the 
disk and the black hole play a lesser role. 



model than the GRFFE model for otherwise similar mod- 
els. This is likely due to a small extra pressure support 
by the corona or coronal wind. Both polar axes show the 
same behavior, so this is not likely due to residual time- 
dependence in the GRMHD model. Despite this effect, the 
opening angle for most of the power output at large dis- 
tances is similar for the GRMHD and GRFFE models. This 
means in the GRMHD models that the coronal wind pro- 
vides some additional collimation, but the primary confine- 
ment of the Poynting-dominated jet is due to the corona. 
In the GRMHD models, the gas+magnetic pressure of the 
corona that contains no ordered field is found to be similar 
to the pressure provided by the force-free field in GRFFE 
models. This agreement between the pressures explains why 
the funnel field geometry is similar, but the process that 
sustains this balance is unknown and should be studied. 



5 GRMHD DISK MODEL WITH INITIAL 
v = 3/4 GRFFE FIELD 

In lMcKinnev fc Naravanl (120061 ) we showed that the v = 3/4 
power-law model represents well the currents measured in 
fully turbulent GRMHD numerical models. We have also in 
this paper shown that the magnetic field geometry and the 
power output in the Poynting-dominated funnel region near 
the poles of the black hole are properly modelled by the 



v = 3/4 rotating force-free model. However, in the GRMHD 
models there are no ordered field lines that thread the disk 
or corona. What happens to these field lines? 

In this section we study the process whereby the 
turbulence-driven convective and magnetic instabilities in 
the disk and cor ona destroy an ordered field. T his process 
was discussed in iMcKinnev fc Gammiel <|2004h ; iMcKinnevI 
(2005a|), where it was indicated that the final geometry of 
the field lines was nearly independent of the initial field ge- 
ometry. The initial field geometries in those studies include 
a single field loop, many field loops with different poloidal 
directions, and a uniform net vertical field. In all cases, the 
final field geometry within the accretion disk, corona, and 
around the black hole was qualitatively similar (see also 
lMcKinnev||2005al ). Here we perform a similar study with the 
non-rotating v = 3/4 solution as the initial field geometry. 



5.1 Numerical Setup 

The GRMHD equations of motion are integrated numeri- 
cally usi ng a modified version of a numerical code called 
HARM (G ammie et a l. 2003a), which uses a conservative, 
shock-capturing scheme. Compared to the original HARM, 
the inversion of conserved quantities to primitive vari- 
ables is performed by using a new fa ster and more ro - 
bust two-dimensional non- linear sol ver (|Noble et al.ll2006l ). 
A parabolic interpolation scheme (jColella 19841 ) is used 
rather than a linear interpolation scheme, and an optimal to- 
tal variational diminishin g (TVD) third order Runge-Kutta 
time stepping (|Shu| [l997) is used rather than the mid-point 
method. For the problems under consideration, the parabolic 
interpolation and third order time stepping method reduce 
the truncation error significantly, even in regions where 
b 2 /p > 1. 



5.2 GRMHD Model with v = 3/4 Initial Field 

The v = 3/4, a/M = solution for the poloidal field 
is implanted in th e same torus model o f the accretion 
disk as studied in IMcKinnev fc Gammiel (|2004l ) and the 
G RMHD equations are integrat ed. As in the fiducial model 
of IMcKinnev fc Gammiel (|2004l ). the black hole has a spin 
of a/M = 0. 9375, which is close to the equilibrium value of 
a/M ~ 0.92 (|Gammie. Shapiro, fc McKinnevlliooi ). 

As a reference point, Figure [TBI shows the initial and 
final field geometry for the GRMHD model with a/M = 
0.9375 with an initial field c orresponding to a single field 
loop in the torus, as studied bv lMcKinnev fc Gammiel (|2004T ) 
and others. 

Figure [T7] shows the initial and final field geometry 
for the GRMHD model with a/M = 0.9375 that uses the 
v = 3/4 force- free solution for the initial field geometry. 
A model initially with a net vertical field leads to qualita- 
tively similar results. The GRMHD model with a single loop 
exhibits more turbulence and larger variations in the polar 
enclosed poloidal current, which at high/low levels saturates 
to red/black in the color figure. 

The final poloidal field distribution at the poles follows 
the v = 3/4 force- free solution for both the models as shown 
in Figures [TBI and [T71 Other models, not shown, such as one 
with an initial net vertical field also agree with the v — 3/4 
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Figure 16. Initial (left panel) and final (right panel) poloidal 
magnetic field geometry (solid black lines) and toroidal field 
strength {Bq, red positive, blue negative) for a GRMHD numeri- 
cal model that was initiated with a gas torus containing a poloidal 
magnetic loop. The model has a rapidly spinning black hole with 
a/M = 0.9375. The panel on the right shows the standard result 
for a GRMHD numerical model of an accretion disk, in which the 
final system has a disk, corona and outflowing wind with disor- 
dered field, and a polar region (the Poynting-dominated jet) with 
highly ordered field. 




jc f/iaw.1 n l b tdm v 

Figure 17. Similar to Figure [TBI except that the initial magnetic 
field corresponds to the v = 3/4 force-free solution. Despite the 
very different initial conditions, the final states in the two models 
are similar. 



force- free solution in the polar region. The value of Qf in 
the disk and transition region is also similar. This is despite 
the fact that the detailed structure of the magnetic geome- 
try in the disk is different. For every model we have studied, 
the turbulence due to the magnetorotational instability effi- 
ciently redistributes the currents to the v = 3/4 power-law 
distribution. 

For these models, the time evolution shows that the 
global magnetic field bifurcates into two regions. The funnel 
region forms as material with low angular momentum falls 
into the black hole and launches a jet, while many of the 
ordered fields threading the torus are tied to material that 
is gaining angular momentum and so is launched backwards 
along the disk surface. 

As the bifurcation unfolds, the pressure support of the 
v — 3/4 field is lost in the region outside the funnel. As the 
disk accretes, some of the disk surface material resupplies 
that evacuated region to replace the pressure that was pro- 
vided by the field. This process forms the corona and drives 
the disk wind. 

A measurement of the gas+magnetic pressure in the 
coronal region shows that the evolved GRMHD model with 
the corona has a similar pressure as the initial force-free 
fields in that region. There is a small additional pressure 
driven by the flux of mass, and this explains why the 
GRMHD model is slightly more collimated than the force- 
free v = 3/4 model. 

A similar quasi-stationary state of the accretion disk 
and field is found with various other initial field geometries 
jHirose et al.ll2004l ; lMcKinnevll2005a! ). Only a contrived ge- 
ometry with a purely toroidal field leads to significant de- 
viations in the quasi-stationary state, and of course models 
with a purely toroidal field have no toroidal currents. 



6 LIMITATIONS 

The primary limitation of the present study is that the 
numerical models are axisymmetric. Therefore, while the 
GRFFE simulations show that the Poynting-dominated jet 
is stable to axisymmetric pinch perturbations, we cannot tell 
whether it is kink stable. 

Another limitation is that in the GRFFE models the 
disk is assumed to be infinitely thin. This is in contrast to 
the GRMHD models where the disk is thick resulting in 
some power from the black hole being absorbed by the disk 
and corona. However, because the force-free region in the 
GRMHD simulations is sharply defined, a good estimate of 
the power can be obtained from the GRFFE models by lim- 
iting the integration of the angular power density to a cone 
around the poles that is the same angular size as that found 
in the GRMHD models. In this way, the Poynting-dominated 
jet power can be estimated for any disk thickness. 



7 CONCLUSIONS 

Motivated by the simple toroidal current distribution we 
found in GRMHD numerical mode ls of turbulent accretion 
disks ( McKinnev fc Naravanl 120061 . in this paper we have 
studied simple force-free (GRFFE) models of the black hole 
and disk magnetosphere. We idealized the GRMHD solution 
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for the accretion disk as an equatorial rotating conductor 
with boundary conditions defined by an equatorial toroidal 
current and field angular frequency. 

To be gin with, in § 3 we conside red the paraboloidal 
problem of iBlandford fc Znaiekl [| 19771 ) in which the current 
is of the form dl^/dr oc r _1 . For a slowly rotating black 
hole and a non-rotating disk, we obtained excellent agree- 
ment between the numerical GRFFE solution and BZ's an- 
alytical solution; this is not surprising since the BZ solution 
was derived in this limit. We then considered a rapidly ro- 
tating black hole with a non-rotating disk and also a rotat- 
ing hole with a Keplerian disk. In these cases, the GRFFE 
model deviates to some extent from the BZ solution, but 
the agreement between the two is still reasonably good. For 
all models, collimation was found to be primarily due to the 
poloidal field geometry determined by the toroidal current 
in the disk. Forces due to hoop-stresses nearly cancel cen- 
trifugal forces. 

In paper I, we considered a self-similar force- free 
model of an equatorial current sheet with a toroidal cur- 
rent following dl^/dr oc r ~ 5 ^ 4 (the v = 3/4 model) 
ijMcKinnev fc Naravanl[2006l 'l. In that paper we showed the 
structure of the field from the force-free solution agrees well 
with the poloidal field in the jet region of GRMHD numer- 
ical models. However, because the force-free model had no 
rotation, it had no toroidal field and we were unable to com- 
pare the toroidal structure or the Lorentz factor of GRMHD 
jets. 

In § 4, we considered a similar force-free model with 
dl,j,/dr oc r _5//4 in the equatorial plane, but with an ar- 
bitrarily rapidly rotating black hole and a Keplerian disk. 
One well-known subtlety of force-free models of the black 
hole and disk is how to define quantities within the tran- 
sition region between the black hole and accretion disk 
ijBlandford fc Znaieldlllmh . We experimented with two dif- 
ferent models, motivated by GRMHD simulations, for the 
field angular velocity in the transition region. 

The agreement between these improved v = 3/4 
GRFFE models and the jet /funnel region of the GRMHD 
models is excellent. We found that the poloidal struc- 
ture of the field in the rotating GRFFE model is nearly 
the same as in the non-r otating model considered in 
iMcKinnev fc Naravanl ((2006) , and like the latter agrees well 
with the poloidal structure of the field in the GRMHD jet. 
Within the jet driven by the black hole, forces due to hoop- 
stresses nearly cancel centrifugal forces, suggesting that self- 
collimation is not important. The toroidal field structure 
in the rotating GRFFE model agrees very well with the 
toroidal structure of the GRMHD jet, as do the power out- 
put in the jet and the Lorentz factor of the jet. Thus, the 
rotating GRFFE model gives an accurate description of the 
acceleration, collimation, jet opening angle, field pitch an- 
gle, electromagnetic power, etc., of the relativistic jet seen 
in GRMHD simulations. 

While there is excellent agreement in the jet regions 
of the two solutions, there are large differences between the 
force- free and GRMHD models in other regions. In the force- 
free model the primary agent responsible for collimating the 
jet is the poloidal field threading the disk. These field lines 
lie external to the field lines in the jet (which thread the 
black hole), and it is the pressure of this external field that 
produces the collimation. In the GRMHD numerical model, 



on the other hand, the region external to the jet is filled 
with a disorganized, magnetized corona and wind, and it is 
the pressure from this gas that causes the collimation. This 
suggests that the corona is required for the collimation of 
the Poynting-dominated jet in GRMHD simulations. 

We also found that the GRFFE jet is stable to axisym- 
metric Eulerian perturbations. In the GRFFE model, the 
black hole jet is in force-balance with the fo r ce-free disk 
wind. In the GRMHD models of IMcKinnev! (|2006d\ we 
found that at large radii the jet is unstable to axisymmet- 
ric perturbations. This is consistent with the lack of coronal 
material (and so coronal pressure support that replaces the 
magnetized support in the force-free regime) at large radii. 
This suggests that astrophysical jets require support by a 
corona or disk wind in order to remain stable, but at large 
radii such support may be naturally unavailable and this 
may lead to toroidal-fi eld-driven instabili ties and dissipation 
in the jet as shown in lMcKinnevi (|2006d ). 

The lack of an ordered field in the corona is related 
to the large difference between the GRFFE and GRMHD 
models with respect to the power output from the disk. 
In the force-free model, the disk field lines carry out elec- 
tromagnetic power in the form of Poynting flux. In the 
GRMHD model, however, little electromagnetic power from 
the disk reaches large radii. Instead, most of the electro- 
magnetic energy from the disk is dissipated in the mat- 
ter and converted to other forms of energy in the plasma 
IMcKinnev fc Naravarj|2006l '). 

In § 6, we described a numerical experiment in which 
we further explored the connection between the v — 3/4 
ordered GRFFE solution and the magneto-rotationally un- 
stable turbulent GRMHD model. An accretion disk was em- 
bedded with the v = 3/4 ordered field and followed using 
the GRMHD code. During the evolution of the system, the 
toroidal current in the disk retains its initial v = 3/4 behav- 
ior. However, the large-scale fields associated with the ini- 
tial state are driven outwards by angular momentum trans- 
port and the pressure of those fields is replaced by coronal 
gas+magnetic pressure. This suggests that there is a con- 
tinuous coupling between the currents in the disk and the 
coronal material. The corona and coronal wind are gener- 
ated by reconnection processes and magnetic stresses in the 
disk. Even without an orga nized field the coronal win d is 
collim ated, as suggested by iHeinz fc Begelmanl (|200fj| ) ; ILl 
(2002). The gas+magnetic pressure in the corona is similar 
to the magnetic field pressure in the GRFFE model, and it is 
because of this that the jet regions are so similar in the two 
models. It remains unknown by what mechanism the pres- 
sure of the corona in GRMHD models drives the power-law 
v — 3/4 toroidal current that generates the corresponding 
magnetized jet. 
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APPENDIX A: ASYMPTOTIC FORM OF 
PITCH ANGLE AND LORENTZ FACTOR 

In this appendix, we show how the Lorentz factor scales with 
radius for a given power-law distribution of currents in the 
disk. First, we remind the reader that the general relativis- 
tic language can be translated into an special relativistic 
language (3+1 splitting) by defining a space-time foliation. 
A convenient foliation is to choose to measure all quanti- 
ties with respect to an observer with 4-velocity of a zero 
angular momentum observer (ZAMO) frame. This is given 
by = {~ a > 0,0,0}, where a = l/y/—g u . In this frame, 
the force- free equations take the same form as in Minkowski 
space-time. 

In the ZAMO frame the drift velocity is written as 

i _ e^ k EjB k 



B- 



(Al) 



where e 8:,fc is the spatial permutation tensor. The Lorentz 
factor is given by 

(A2) 



as in special relativity. The electric field is given by 

E l = -t^VjBk, (A3) 



Given flux-freezing, axisymmetry, and stationarity, the fluid 
velocity or drift velocity is related to the magnetic field and 
field angular frequency via the relation 

Rfl F 



B„ 



B t 



(A4) 



where these are orthonormal basis quantities where p stands 
for p oloidal, t stands for toroi dal, and R is cylindrical ra- 
dius (Bck enstein fe O ron 197j j). For a review of the force- 
free equations see McKinnevI ( 2006aT ). 

Now, consider the asymptotic regime where r — > oo, 
and assume that the toroidal field dominates in this regime 
such that the mostly radial poloidal velocity approaches the 
speed of light. Then from equation ()A3[) one finds that 



\E\ 



\E 



\B t 



(A5) 



From the drift velocity, one has that 



2 IE x B 

V = 



E' 



B\Bl + BlY (B% + B1Y (Bp/St) 2 + L 
and so the Lorentz factor can be written as 

B v 



(A6) 



(A7) 



Beyond the light cylinder, v t <C RQ,p and assuming v p ap- 
proaches the speed of light, equation (|A4[) then gives that 

|i = -R£l F . (A8) 
Bp 

Notice that this result is consistent with the numerical re- 
sults of this paper, where we discuss the pitch angle such as 
given by equation (|15p . 

The above implies that 



T = RSlp = Q F r 



(A9) 



for a field geometry that follows 9j ■ oc r a . In general, for 
a solution with a current distribution with power-law index 



v, in the slow-rotation approximation, the minimal-torque 
solution (|MichelH l969) gives an opening angle of approxi- 
mately 



and so 



r oc r 



-v/1 



l-v/2 



(A10) 
(All) 



For example, for the paraboloidal field geometry with v 
1, one has that o 



-1/2 



= 3/4 one has that 9j 



so that F oc 
375 so that T 



r 1/2 . For 

„ „0.625 



This basic scaling is fairly consistent with the results of 
this paper, where deviations are expected for a general 
black h ole spin and d isk ro t ation profile. See also |Beskin 



d 19971); iBeskin et 
iBeskin fe Nokhrinal 



ail (Il998| 
rd20Q6|);t 



I); IBeskin fc Malvshkinl (|2000h 



Naravan et all (|2006l ) 
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